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Experimental  observations  of  the  time  and  space  evolution  of  resonantly 
enhanced  electric  fields  and  plasma  density  in  cylindrical  geometry  demon¬ 
strate  the  development  of  caviton  structure  in  the  direction  perpendicular  to 
the  driving  electric  field.  Electrostatic  fields  in  this  perpendicular 
direction  are  observed  to  have  growth  times  on  the  ion  time  scale  and  to 
develop  concurrently  with  two-dimensional  density  profile  modification  near 
the  critical  surface. 
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SIMULATION  OF  IONOSPHERIC  EM-PLASMA  INTERACTION  IN  A  LARGE  LABORATORY  DEVICE 
(Development  of  2-D  Cavitons)  "  ~  Final 

Introduction 

The  interaction  of  electromagnetic  (EM)  waves  with  plasma  is  an  area  of 
research  relevant  to  such  diverse  applications  as  laser-pellet  fusion,  RF/ 
microwave  heating  of  plasmas,  and  HF-induced  ionospheric  modification.  In 
these  applications  one  is  interested  in  transforming  the  incident  EM  radiation 
into  particle  energy,  thus  heating  the  plasma.  In  a  collisionless,  unmaanetized, 

g 

inhomogeneous  plasma  this  transformation  occurs  in  the  following  way  .  Propa¬ 
gating  EM  waves  with  frequency  w  are  reflected  at  the  point  where  n  =  ncrit  Cos2e 
where  ncri-t  is  the  density  where  ui  =  o)pe.  At  higher  densities  the  wave  is 
evanescent  and  the  amplitude  decays  exponantially  in  space.  For  p-polarized 
EM  waves  (see  Fig.  1)  the  component  of  E  parallel  to  the  density  gradient  can 

drive  an  electrostatic  resonance  at  n _ .. .  This  electrostatic  resonance  gives 

crit 

rise  to  electron  plasma  waves  which  are  then  Landau  damped,  thus  heating  the 
electrons. 

Since  this  process  occurs  near  the  critical  density  EM  propagation  effects 
in  the  subcritical  plasma  do  not  contribute  to  the  physics  of  absorption.  One 
can  therefore  simplify  the  experiment  further  by  modeling  the  evanescent  elec¬ 
tric  field  between  EM  cutoff  and  the  critical  point  with  an  electrostatic  RF 

1  3 

field  —  the  so-called  capacitor  plate  type  of  experiment  ’  . 

Up  until  now,  most  microwave  and  capacitor  plate  investigations  of  critical 
surface  phenomena  have  dealt  only  with  phenomena  in  the  dimension  defined  by 
the  direction  of  the  initial  density  gradient;  the  plasma  density  in  the  direc¬ 
tions  perpendicular  to  the  density  gradient  was  assumed  to  be  uniform  both 

initially  and  as  the  profile  evolved.  However,  computer  simulations  (both 
7  8 

fluid  and  kinetic)  and  laser-plasma  experiments  have  indicated  that  two-  and 
three-dimensional  effects  can  play  an  important  role  in  EM-plasma  interactions. 
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Furthermore,  2-D  numerical  calculations  on  the  evolution  of  Langmuir  waves  have 

5  6 

shown  that  a  1-D  soliton  can,  if  perturbed,  collapse  in  two  and  three  dimensions  '  . 

Summary  of  work  performed 

Stimulated  by  the  above-mentioned  results,  we  have  performed  experiments 
designed  to  investigate  the  two-dimensional  evolution  of  an  electrostatically 
driven  inhomogeneous  plasma.  In  pursuing  this  work  we  have 

Developed  a  new  cylindrically  symmetric  experimental  device  especially 
suited  for  a  2-D  investigation. 

Developed  2-D  diagnostics  which  allow  us  to  measure  phenomena  anywhere 
on  the  2-D  plane. 

Set  up  a  computerized  data  taking  system  to  efficiently  handle  the 
great  volume  of  data  involved  in  a  2-D  investigation. 

Made  the  first  observation  of  the  development  of  two-dimensional 
structure  in  cavitons.  These  observations,  which  have  been  submitted 
for  publication,  are  described  in  detail  in  the  next  section. 

2 

Developed  a  simple  theoretical  model  based  on  previous  1-D  work  to 
explain  our  results. 
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Development  of  Two-Dimensional  Structure  in  Cavitons 

The  experiment  was  performed  in  the  cylindrically  symmetric  device  shown 

in  Figure  2.  A  weakly  ionized  Argon  plasma  (ng  =  109  cm’3,  =  2  x  10’1*  torr) 

was  produced  by  a  pulsed  HF  oscillator  (v  =  10  MHz,  pulse  length  *  25  vs,  duty 

cycle  <  25%)  connected  to  a  long,  cylindrical  antenna  consisting  of  twelve  rods 

extending  axially  and  equally  spaced  azimuthally.  The  experiment  is  performed 

in  the  quiescent  discharge  afterglow  where  <5n/n  <  0.3%,  Tg  =  1.0  eV,  T..  =  0.1  eV, 
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and  Xp  =  0.23  mm.  The  density  gradient  in  the  experimental  region  is  primarily 

1  Zn  \ 

in  the  radial  direction  [L  =  —  —  1" 1  =  16  cm]  and  the  coaxial  line  formed 

,  no  dr  / 

by  the  central  rod  and  the  chamber  wall  provides  a  purely  radial  quasistatic 
driving  field  (vQ  =  285  MHz,  P  <  100  watts,  Eyac  <  4  V/ cm  typically)  to  produce 
the  caviton.  This  is  in  contrast  to  the  1-D  capacitor  plate  experiments^  where 
fringing  driver  fields  in  the  direction  perpendicular  to  the  density  gradient 
are  unavoidable.  The  cylindrical  geometry  also  allows  one  of  the  directions 
perpendicular  to  the  density  gradient  to  close  on  itself,  thus  avoiding  boundary 


effects . 

The  caviton  electric  field  is  measured  with  an  unperturbing  diagnostic 
electron  beam*  (<J>b  =  7.5  -  9.0  kV,  1^  =  0.1  uA)  with  a  temporal  (spatial) 
resolution  of  0.2  us  (1  inn) .  The  beam  can  be  mechanically  moved  both  radially 
and  azimuthally  which  allows  us  to  diagnose  the  entire  r,  0  -  plane.  The 
plasma  density  is  measured  with  a  small  wire  probe  (10  mil.  diameter,  0.2” 
long)  positioned  by  a  two-dimensional  probe  drive.  Movement  of  the  probe  is 
controlled  by  a  computer  which  also  stores  the  data  on  flexible  magnetic  disks. 
Time  resolution  is  obtained  through  the  use  of  a  boxcar  averager. 

The  basic  experimental  sequence  is  as  follows:  at  a  selected  time  in  the 
discharge  afterglow  a  short  (1  -  10  us)  VHF  burst  (P  <  100  watts)  is  applied 
to  the  center  electrode  of  the  chamber  to  drive  the  caviton.  At  a  selected 
time  t  within  the  VHF  burst  a  shorter  (<  1  us)  voltaqe  pulse  is  applied  to  the 
electron  gun.  The  cavi ton-modulated  image  of  this  electron  burst  is  then  mea¬ 
sured  on  the  phosphorus  screens  at  a  selected  coordinate  (r,  8).  Only  fields 
greater  than  4  V/cm  produce  an  observable  modulation  thus  making  this  diag¬ 
nostic  sensitive  to  resonantly-enhanced  fields  (E  ~  100  V/cm)  but  not  to  the 
driver  field  (E  <  4  V/cm).  After  measuring  the  axial  dependence  of  the  caviton 
one  can  unfold  this  image  data  to  get  the  amplitude,  direction,  and  relative 
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phase  of  the  enhanced  VHP  electric  field  at  (r,  8,.t).  When  density  data  are 
taken,  the  probe  current  is  sampled  approximately  2  us  after  the  end  of  the  VHF 
burst  to  avoid  spurious  probe  effects. 

In  previous  experiments^  the  caviton  electric  field  was  always  in  the 
same  direction  as  the  driving  field.  Such  a  response  is  shown  in  Fig.  2b 
where  the  caviton  field  is  purely  radial  and  has  a  shape  characterized  by  a 
large  resonance  and  several  smaller  subcritical  peaks;  this  shape  is  quite 
similar  to  the  Airy  function  solution  of  the  so-called  capacitor  plate  pro¬ 
blem2.  The  plasma  density  response  to  the  ponderomotive  force  of  the  field 
is  shown  in  Fig.  2c.  Our  main  new  finding  is  that  this  caviton  can  evolve 
in  such  a  way  as  to  produce  a  strong  electric  field  component  in  the  direc¬ 
tion  perpendicular  to  the  driving  field. 

Figures  3a,  b,  c,  show  the  magnitude,  direction,  and  relative  phase  of 
the  electric  field  at  time  t  *  0.5  us  *  3.3  1  for  a  section  of  the  polar 

plane.  -  Note  that  at  this  early  time  in  the  VHF  burst  the  enhanced  electric 
field  is  purely  radial  and  that  it  has  very  little  azimuthal  dependence.  Two- 
dimensional  density  measurements  show  that  the  initial  (zeroth-order)  density 
has  good  symmetry  and  has  no  strong  gradients  in  the  azimuthal  direction.  Thus 
initially  we  have  a  fairly  symmetric  caviton  that  forms  a  ring  concentric  with 
the  center  conductor  and  has  purely  radial  electric  fields. 

Figures  3d,  e,  f  show  the  same  polar  section  at  a  later  time  in  the  VHF 
burst  (t  =  6.75  us  =  44.6  uip..'1).  The  electric  field  magnitude  has  now  devel¬ 
oped  an  azimuthal  dependence  thus  destroying  the  initial  azimuthal  symmetry. 

The  field  has  also  developed  a  vector  component  in  the  ©-direction  that  is  com¬ 
parable  in  maanitude  to  the  radial  component.  More  detailed  temporal  measure¬ 
ments  show  that  the  growth  and  decay  time  for  |Ej  is  -  1.5  us  =  10  u_._1.  As 

o  pi 

indicated  by  the  Lissajous  patterns  in  Fig.  3e ,  there  is  a  definite,  spatially 
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varying  phase  relationship  between  the  E  and  EQ  components  of  the  measured 

r  v 

electric  field. 

Two-dimensional  plasma  density  measurements  are  shown  in  Figure  4.  Ex¬ 
perimental  conditions  are  similar  to  those  existing  during  electric  field  measure 
ments.  The  critical  density  occurs  at  r  =  8  cm  and  consecutive  contours  cor¬ 
respond  to  a  density  change  of  — =  1.26%.  Note  that  a  larger  polar  section 

ncrit 

is  shown  here  than  in  Figure  3.  Figure  4a  shows  the  zeroth  order  (undriven) 

density  which  decreases  with  increasing  r  and  is  fairly  symmetric  azimuthally. 

The  density  ripples  at  large  r  are  due  to  the  regularly  spaced  rods  of  the  plasma 

production  antenna.  At  the  critical  surface  the  magnitude  of  these  ripples  is 

—  =  0.5%  with  an  azimuthal  wavelength  of  X  =  4.1  cm  =  180  X„. 
n  D 

Figures  4b  and  4c  show  isodensity  contour  plots  taken  after  a  1  ys  and 
8  ys  VHF  burst,  respectively.  Dotted  contours  indicate  a  closed  contour 
which  contains  a  local  minimum  (density  hole).  As  with  the  electric  field 
measurements,  we  again  see  an  increase  in  the  azimuthal  asymmetry  with  VHF 
burst  length.  In  Figure  4c  radially  and  azimuthally  localized  density  struc¬ 
tures  have  developed.  Note  that  the  azimuthal  spacing  of  these  structures 
(A6  =  30°)  matches  the  azimuthal  peak-to-peak  spacing  of  the  enhanced  electric 
field  (Fig.  3d,  f).  Further,  note  that  the  azimuthal  locations  of  these  density 
structures  correspond  to  the  density  minima  of  the  perturbation  caused  by  the 
rods  of  the  plasma  production  antenna. 

A  simple  model  to  explain  the  development  of  two-dimensional  structure 
from  an  initial  perturbation  is  as  follows.  Consider  a  linear  density  profile 
with  a  periodic  perturbation  in  botn  the  density  and  the  density  gradient  scale 
length  L  =  no(3n/3x)_1: 

n(x,v)  =  nQ  (1  +  x/l(y)  +  eCosky)  (1) 

where  nQ  is  the  critical  density  and  L(y)  =  LQ  (1  +  eCosky)”1.  Here  x  is  the 
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direction  of  the  drivinq  electric  field  (correspondino  to  our  r)  and  y  is  perpen- 

2 

dicular  to  the  driver  (our  0).  Adopting  the  scaling  of  Morales  and  Lee  ,  the 
solution  for  E  when  3E  /3x  »  3E  /3y  is  E  (x,y)  =  (k-1)2/3  E  A(n,  t)  where 

A  A  A  A  U  « 

n  =  (kpL)2/3  (x/L  +  eCosky),  kQ  =  (Debye  length)-1,  EQ  is  the  driver  field,  and 

A  is  the  spatial  variation  of  the  enhanced  electric  field  and  has  amplitude  of 

order  one.  Using  the  electrostatic  condition,  3E  /3y  =  3E  /3x,  and  integrating 

x  y 

over  x  one  can  solve  for  Ey.  Near  the  critical  surface  and  to  first  order  in  c 
we  obtain: 

Ey(x,y)  =  -  (Lq  ek  Sinky)  Ex(x,y)  (2) 

The  critical  surface  with  nonlinear  perturbation  is  now  defined  by  x/L  +  eCosky 

+  6n/n  *  0  where  <5n/n  =  -  E2/16imT  and  E2  =  E2  +  E2.  We  can  rewrite  this 
o  o  0  6  x  y 

as  n  =  Y  P(y)  A2(n,t)  where  y  =  (koL0E0)  Z/16lTn0Te>  p(y)  =  (1  +  5  Sin2ky)/ 

(1  +  eCoskv)2,  and  6  =  (L^k)2.  We  must,  in  general,  solve  this  transcendental 
equation  for  n  for  each  value  of  y.  However,  since  our  perturbation  is  small 
P(y)  a  1  and  we  can  expand  A2(n,t)  around  the  value  nQ  =  Y  A2(no,t).  Keeping 
terms  of  first  order  in  e  and  6  we  obtain  for  the  modified  profile  near  the 
critical  surface: 

n(x,y)  *  nQ  ^  1  -  p  +  x/L(y)  +  [1  +  y(4/3  +  28)]  eCosky 

-  y(l  +  6)  6  Sin2k.y^  (3) 

where  8  =  2nQ  A' (n0,t)/[A(n0,t)  -  2nQ  A ' (n0t ) 3  and  y  =  n0/(kDL0)2/3.  If  we 

choose  for  A(n,t)  the  steady  state  Airy  function  solution  it  can  be  shown  that 
for  all  positive  n0«  -0.2  £  B  <  0.  Then  comparing  eouations  (1)  and  (3)  one 
can  see  that  the  nonlinear  modification  to  the  density  near  the  critical  sur¬ 
face  has  increased  the  amplitude  of  the  density  perturbation.  In  addition, 
second  harmonics  of  the  initial  perturbation  are  produced  by  the  Sin2ky  term 


which  represents  the  beginning  of  localization  in  the  y-direction.  Thus  this 
approximate  model  suggests  a  mechanism  for  the  production  of  an  E^  electric  field 
component  with  a  spatially  varying  phase  relationship  to  Ex  and  for  the  non¬ 
linear  development  of  two-dimensional  density  structure  from  an  initial  pertur¬ 
bation. 

5 

Our  observations  are  also  reminiscent  of  theoretical  studies  of  soliton 
collapse  in  two-dimensions.  Onaoing  theoretical  work^  has  shown  that  a  similar 
collapse  can  occur  in  an  inhomogeneous,  externally  driven  plasma.  More  de¬ 
tailed  measurements  of  the  structure  growth  rates  and  the  traoping  of  VHF  fields 
should  reveal  which  model (s)  should  be  used  to  explain  our  data.  Our  work  may 
also  be  relevant  in  explaining  the  small  scale  rippling  of  the  critical  surface 
as  observed  in  laser-plasma  simulations  and  experiments7,  although  such  rippling 
is  usually  attributed  to  electromagnetic  effects. 
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Figure  1 . 


Figure  2a. 


Figure  2b. 


Figure  2c. 
Figure  3a. 


Figure  3b. 

Fiqure  3c. 

Figure  3d. 


Fiqure  3e. 


Figure  Captions 

Schematic  drawinq  of  the  propagation  of  a  p-polarized  electro¬ 
magnetic  wave  in  an  inhomogeneous  plasma. 

Schematic  diagram  of  the  experimental  apparatus  showing  coaxial 
experimental  area  on  right  and  movable  diagnostic  beam  on  left. 
Caviton  field  amplitude  vs.  radial  position  at  t  =  0.5  ys  and 
azimuthal  position  9  =  90°  (0°  being  top  of  chamber  and  increasing 
clockwise  facing  phosphor  windows).  Enhanced  fields  at  inner 
conjugate  resonance  point  (r  =  3  cm)  are  outside  the  range  of  the 
electron  beam  diagnostic. 

Electron  density  profile  corresponding  to  field  profile  in  Fig.  lb. 
Contour  plot  of  caviton  E-field  amplitude  at  t  =  0.5  vs  for  a  section 
of  the  polar  plane.  Contour  level  is  proportional  to  field  ampli-. 
tude  and  level  changes  by  6  V/cm  with  each  line.  Minimum  and 
maximum  values  are  zero  and  144  V/cm  respectively. 

Equal  length  arrow  plot  for  t  =  0.5  ys.  Direction  of  arrow  gives 
direction  of  electric  field. 

Three  dimension  graph  of  the  same  data  as  in  2a.  Third  dimension 
is  proportional  to  |Ej. 

Contour  plot  of  total  caviton  electric  field  at  t  =  6.75  ys  for 
same  polar  section.  Contour  level  changes  by  3  V/cm  with  each 
line.  Minimum  and  maximum  values  are  zero  and  85.5  V/cm  respec¬ 
tively. 

Equal  length  arrow  plot  for  t  =  6.75  ys.  Relative  size  of  Er 
and  Eq  can  be  determined  from  inclination  of  arrow  and  relative 
phase  can  be  determined  from  arrow  direction.  When  relative 
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phase  was  90°  an  ellipse  was  drawn  with  (major  axis)/(minor  axis) 

|Erl/|E0l- 

Figure  3f.  Three  dimensional  graph  of  the  same  data  as  in  2d. 

Figure  4.  Isodensity  contour  plots  a)  before  the  VHF  burst,  b)  after  a 
1  us  VHF  burst,  and  3)  after  a  8  vs  VHF  burst. 
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The  following  is  the  title  and  abstract  of  our  work  performed  under 

which  was  presented  at  the  American  Physical  Society 
Meeting,  San  Diego,  California,  November  10-14,  1980. 


9T  20  Development  of  Two- Dimensional  Structure  in  Cav- 
]  tons.*  D.  1.  EGGLESTON .  K.  T.  »ONC.  C.  »■  DARROW.  UCLA. 
--A  coaxial  capacitor  plate  experiment  in  which  tha  dn* 
ving  field  and  the  density  gradient  are  both  in  the  rad¬ 
ial  direction  only  is  performed  to  study  RF-plasaa  inter¬ 
actions  in  two  dimensions  (v  ■  28S  fftz).  Tha  time  and 
space  evolution  of  resonantly  enhanced  electric  fields 
and  plasma  density  are  measured  in  two-dimensions  by  a 
non-perturbing  scanning  diagnostic  electron  beam  and  a 
thin  wire  probe,  respectively.  These  measurements  dem¬ 
onstrate  the  development  of  caviton  structure  in  the  di¬ 
rection  perpendicular  to  the  driving  electric  field. 
Electrostatic  fields  in  this  perpendicular  direction  are 
observed  to  develop  concurrently  with  density  profile 
modification  and  critical  surface  rippling  and  have 
growth  times  on  the  ion  time  scale.  Possible  mechanisms 
for  this  two-dimensional  evolution  are  considered. 
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